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Abstract 

It is shown that in the large N c limit heavy baryon masses can be estimated quantitatively in a 1 /N c expansion using the Hartree 
approximation. The results are compared with available lattice calculations for different values of the ratio between the square root 
of the string tension and the heavy quark mass yfcr/niQ. These estimates implement important 1 /N c corrections and assume a string 
tension independent of N c . Using a potential adjusted to agree with the one obtained in lattice QCD, a variational analysis of the 
ground state spin averaged baryon mass is performed using Gaussian Hartree wave functions. Relativistic corrections through the 
quark kinetic energy are included. The results provide good estimates for the first sub-leading in 1 /N c corrections. 


1. Introduction 

QCD in the large N c limit becomes a non-trivial theory in 
terms of an arbitrary and fixed t’Hooft coupling A = a s N c (TJ. 
In that limit, baryons j2), unlike mesons, remain as complicated 
structures (for a recent review see e.g. 00 and references 
therein). This is the result of the strong coupling of mesons 
to baryons 0( \/N c ), giving baryons a light meson cloud which 
contributes to its mass at leading order in N c . In the world of 
QCD with only heavy quarks, the meson cloud becomes sup¬ 
pressed in A QCD/mQ. ln Q being the heavy quark mass, and 
baryonic states become amenable to a treatment based on non- 
relativistic QCD. Thus, heavy baryons are a good laboratory to 
study the l/N c expansion. This simpler setting of QCD per¬ 
mits a straightforward application of the mean field approach, 
which will be used in the present work and which should pro¬ 
vide a good description of baryons in the large N c and large 
quark mass limits. 

The quantitative understanding of the l/N c expansion has 
become possible in the light meson sector 0, where meson 
masses have been determined in lattice QCD (LQCD) calcu¬ 
lations at different values of N c and in the quenched approxi¬ 
mation, where the leading 0( 1 /N c ) corrections are absent, and 
moderate N c values allow for a safe extrapolation to the large 
N c limit. In addition, estimates based on short distance con¬ 
straints provide an analytical understanding of those results 0. 
More recently, LQCD calculations of low lying baryon masses 
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for N c = 3, 5 and 7 00 have opened the door for a quantita¬ 
tive test of the 1 /N c expansion in baryons as well. Those pio¬ 
neering calculations, which are in the quenched approximation, 
have quark masses in the light to moderately heavy range. The 
present work is largely motivated by the possibility that such 
LQCD calculations could be extended to heavier quark masses, 
where the framework presented here would become realistically 
applicable. 

In his seminal paper, Witten (2) discussed specifically heavy 
baryons in the large N c limit and invoked the mean field Hartree 
approximation. For heavy quarks, it is built from the simple 
two-body Hamiltonian, where the interaction is the OGE (one 
gluon exchange) (see 0 for details) for the short range part of 
the interaction. In addition, there are the long range confining 
forces, whose effects become suppressed as ihq grows, and also 
short distance radiative corrections must be taken into account 
(running of orQ (see ni). Furthermore, the effects of three- 
body interactions are of potential interest; for a recent discus¬ 
sion in the quark model see Ref. fTTl . They will be discussed 
briefly in this work. 

At leading order in the 1 /N c expansion, the ground state of 
the heavy baryon will be described by a wave function which is 
the direct product of single-quark wave functions. Since the hy- 
perfine interactions have spin-flavor non-singlet effects which 
are O (1 /N c ), it is clear that at leading order the spin-flavor state 
of the ground state baryons is in the totally symmetric spin- 
flavor state, and the baryon has a spin-flavor contracted sym¬ 
metry naira, which holds in the limit N c — » oo at fixed quark 
mass. The effect of removing the center of mass (CM) motion 
is sub-leading in 1 /N c , and can be implemented using standard 
techniques such as the Peierls-Yoccoz projection (for a review 
see e.g. mna and references therein). 


Preprint submitted to Physics Letters B 


July 17, 2015 





The mean field for heavy quarks at large N c has been studied 
in Refs, naira along with possible implications for baryonic 
matter. This work builds on that one and compare to recent lat¬ 
tice calculations for /V, = 3, 5, 7 00 after including some 
important 1 /N c effects such as the CM correction. Brief discus¬ 
sions of the role of hyperfine splittings as well as the expected 
corrections of many-body forces are also given. A previous 
large N c analysis has been conducted in Ref. Ifl8l . 

Note that in order to have low lying baryons with different 
spins it is necessary to have more than one flavor of heavy 
quark. The mass of the baryon will then have an O (1 /N c ) hy¬ 
perfine contribution (dependent on the spin S of the baryon). 
The masses of ground state baryons take the form of a rota¬ 
tional band, 

M b (S) = N c m o + ^(S(S + 1) - jN c ) + 0(l/N>), (1) 

where mo and Chf are 0(n®\ and have an expansion in 1 /N c , 
and depend on the quark mass mg. The hyperfine independent 
component of the baryon mass given by mo is obtained by the 
following combination of baryon masses: 

2 

17lc\ — - 

N 2 (N C + 1 )(N C + 3) 2 

Nc_ 

2 

x J^(3+N c (3N c + 2)-m c -3)S)M B (S). (2) 
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The baryon masses studied here will be the ones with the 
hyperfine effects removed, i.e.. Mg = N c mo- These will be later 
compared with the available LQCD results of Refs. iTTl f8l ITSlI . 

Of course, for any different value of N c one has a different 
theory. Thus, in order to relate them one must assume that some 
observables are N c independent. Actually, on general grounds 
one has that: 


mo 




= F(N c 



(3) 


where cr sets the scale of QCD and can be identified for instance 
with the string tension, and mg is the heavy quark mass. F is a 
universal function O (/V ( °j which admits an expansion in 1 /N c , 
and which for large nig can be more conveniently expressed as 

The present work goes beyond Refs. inana by analyzing 
the main 1 /N c contributions such as the CM effect, and rela¬ 
tivistic corrections, and actually compares to available LQCD 
results. For N c = 3, triply heavy baryons have been studied on 
the lattice as a Q-bbb state m, and also re-addressed in quark 
models within several schemes mmm which, however, 
have not addressed larger N c values. 

One important goal on the lattice has been to make the quarks 
as light as possible. Actually, quarkonium studies based LQCD 
proceed always through the determination of the QQ potential, 
and a subsequent solution of the non-relativistic Schrodinger 
equation (see e.g. lf2T1 ). The present work takes a similar point 
of view as a A r -hodv problem. It should be emphasized that 
studying heavy baryons at varying values of N c will help with 


the understanding of the 1 /N c expansion in a setting where an 
analytic approach with small model dependencies can be ap¬ 
plied. 


2. Color singlet states 


The starting point is the Hamiltonian for heavy quarks. Using 
non-relativistic heavy quark field operators Q(x), the Hamilto¬ 
nian is given by: 
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-~^—Q\x) A Q(x) + m Q Q f (x)Q(x) 

2m n 


(4) 


+ A | ,fxcPx' QHx)^-Q(x) 2 f (x')y Q(x / )V(x - x '), 


where A a are the S U(N r ) generators in the fundamental repre¬ 
sentation, and in perturbation theory V(r ) = a s /r is the OGE 
interaction. Here, only two-body interactions are included. The 
role of many body interactions is commented below. An equiv¬ 
alent representation for the case of a heavy baryon is the Hamil¬ 
tonian 
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*a( 0 ® A a (j)V(Xj - Xj) (5) 
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The A ® A interaction implies exact Casimir scaling of the po¬ 
tential energy. Casimir scaling for the QQ potential holds per- 
turbatively up to two loops (there are three-loop violations) m 
and numerically on the lattice ll23l . 

For a colour singlet state the wave function is completely 
symmetric in the orbital and spin-flavour quantum numbers, 
and the baryon behaves effectively as a bosonic system. In 
particular, for ground state baryons the wave function is the 
product of a symmetric spacial wave function and a symmet¬ 
ric spin-flavor wave function and reads as follows: 


'T(xi, ...,x N ) = ..., x n )xsf , 


( 6 ) 


where xs f is the spin-flavor wave function. For excited baryon 
states, spin-flavor and spatial mixed symmetry states also oc¬ 
cur. The color matrix elements for arbitrary N c in the ground 
state can be computed as follows. Starting with the quadratic 
Casimir operator for the fundamental representation given by 
(. F a = A a l2) 


Fq ' F q - Fq ■ Fq - 


N;- l 

2 Nc ’ 


(7) 


for a baryon (colour singlet) state one obtains: 


N c 

0 = (B\(YFd 2 \B) 

i= 1 
N c 

= (B\ J(F,) 2 |B) + 2 ■ Fj\B ) 

i=l i<j 

= N c (B\(Fq) 2 \B) + N C (N C - l)(B\Fq ■ P q ,\B), (8) 


and likewise for a meson state one obtains: 


0 = (M\(Fq + Fqf\M) 

= 2(M\(Fqf\M) + 2{M\F q ■ F q \M) (9) 
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These equations lead to 


(B\F q ■ F q '\B) = -Ml+ 3-1 (10) 

(M\P q ■ Pq\M) = —%jf~ dD 

At very short distances the potential between a heavy quark and 
antiquark should be described with perturbative QCD, and ap¬ 
proximately given by an A,.-independent expression at leading 
order (LO) in terms of the running strong coupling a^ c (r). 


V q\ f W - - 


N 2 C - 1 a N fi{r) 1 


2 Nr 


r 11 log(rA^g) 


( 12 ) 


At long distances it is of linear confining form and the corre¬ 
sponding string tension cr is determined in LQCD. For N c = 3 
the QQ potential has been computed in LQCD in the quenched 
approximation ll24l . and for N c > 3 also 018]. For N c = 3, it is 
well described by the bosonic string model ll25l . namely: 


y eb =3(r) = ~Wr +ar ' 


(13) 


The Coulomb term on the RHS is what results from the fluc¬ 
tuations of the string. It is remarkable that it provides the 
bulk of the Coulomb interaction down to the lattice spac- 
ings used in present day calculations. Using A^g-/ yfer = 
0.503(2)(40) + 0.33(3)(3)/A 2 + 0(N; 4 ) obtained in [25] one 
gets that at r \[cr ~ 0.2 the 1/r term in Eqs. (12 1 and (Tsj ) co¬ 
incide. For the heavy quark mass corresponding to Compton 
wave lengths much smaller than present lattice spacings, where 
the long distance potential plays a minor role, the Coulomb in¬ 
teraction will increasingly become the one predicted by pertur¬ 
bative QCD, Eq. ( [T2| . 

At arbitrary N c , V will only receive corrections 0[\/N 2 ^, 
as required by the l/N c expansion in pure gluodynamics. As¬ 
suming the leading scaling in N c for a s and cr, and Eq ( fTT) , the 
potential becomes: 
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V N -\r) 


= (l+tf(l/A 2 ))Vj 3 (r) 


(14) 


This N c dependence will be loosely named ’’Casimir scaling”. 
This is verified by the t’Hooft coupling A — AnN c a s used in 
Refs. 00. Clearly this follows only if the above assumption 
is made, and with the present calculation at N c > 3 it can be 
verified, as discussed below. 

As mentioned earlier, the 1/N C expansion requires defini¬ 
tion because it compares different theories. The most obvious 
way to proceed is to require that certain quantities are indepen¬ 
dent of N c , e.g., the string tension and quark masses at a given 
scale. Since the LQCD results of Ref. 0 0 have the property 
that the string tension is approximately independent of N c , i.e., 
cr = | ^2 ~ct( 3) ~ const, this condition is adopted in what fol¬ 
lows. The result from Fig.[T|vividly shows the N c independence 
of the QQ potential within the current lattice uncertainties and 
the astonishing agreement with the bosonic string model l25l . 


Thus, generalizing the N c lattice findings l24l the potential for 
all N c will be taken to be: 


V Nc -(r ) = V Nc - 2 (r) 
QQ y ’ QQ w 


n 

\2~r 


+ err 


(15) 


From Eqs. (10-15) the two-body interaction potential in the 
baryon becomes: 


V QQ M = 


N n - 1 


1 in \ 

= --1- err 

N c - 1 V 12r / 


(16) 



Figure 1: Quark-antiquark Potential on the lattice in units of the string tension 
for N c = 3,5,1 compared with the bosonic string model (25) (full line). The 
values for different N c '. 3 (blue), 5 (red) and 7 (black), have been transported to 
avoid cluttering of points. 


3. Mean field approximation and beyond 

3.1. Mean field approximation 

The calculation for different values of N c = 3, 5, 7,... of 
the baryon mass with the Hamiltonian Eq. requires solving 
separate few body problems with their inherent technical com¬ 
plications. In the large N c limit, however, an important sim¬ 
plification arises as a mean field approach becomes valid. The 
structure of the spacial wave function is of Hartree type 0. 
and in the particular case of the ground state it reads: 

N 

if/(xu ... ,x N ) = ]~[ <p(xi) . (17) 

i=i 

For a single baryon, the baryon mass Mb = = (FI is 

given by: 

Mb = N c m Q +N c f d 3 x—^— |V0(x)| 2 
J 2 m Q 

+ —— J'd 3 xd 3 x'\(f>(x')l 2 \<p(x)l 2 V Qe (x- x'). 

(18) 

The large N c scaling becomes obvious after the relation, 
Eq. (pA] is used. It is useful to define the effective mean field 
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potential V(x) generated by N r - 1 quarks 

V(x) = (N c - 1) J d 3 x'V QQ (x - x')\(/)(x')\ 2 

= f d 3 x'V Q Q(x-x')\<p(x?) | 2 , (19) 

where the Casimir scaling assumption provided by Eq. ( fl6| ) has 
been used. The mean field potential is the self-energy of a 
quark within the hadron which sees the remaining N c - 1 quarks 
(which are coupled into the anti-fundamental representation F ). 

The mean field equations are then obtained by minimizing 
with respect to a normalized <p(x) leading to the eigenvalue 
problem: 

- —V 2 0(x) + V(x)f(x) = ef(x). (20) 

2m q 


3.2. Numerical and variational solution 


The mean field equations Eqs. (20 1 and (19 1 can be solved 
by iterations until self-consistency solution is obtained. Actu¬ 
ally, for the case cr - 0 the system can be written as a cou¬ 
pled Schrodinger-Newton equation, which was already solved 
in Ref. l27l . A Gaussian ansatz of the form 


* r) = U* 


-P/b 2 


( 21 ) 


yields a good approximation to this solution and allows for a 
simple analytical discussion. 0 


3.3. CM corrections and mass formula 

One standard and well documented problem of the mean field 
approximation in nuclear physics is the violation of Galilean 
invariance mm which is a symmetry of the starting Hamil¬ 
tonian, Eq. <0 namely the invariance under the boost opera¬ 
tion with velocity v, ,... , xn) —* e lmQV ' 22iX ' i ¥(xi ,..., x,y), 
which implies the energy of the moving system to be given by 
E(P) — M + P 1 IIihqNc where the rest mass differs from the 
inertial mass M + N c hiq. 

Since the interest here is to include 1 /N c corrections in the 
calculation, it is important to build a wave function that is an 
eigenfunction of the momentum. This is achieved by imple¬ 
menting, e.g., the Peierls-Yoccoz projection method [15, 1410 
However, for the simple Gaussian single particle wave function, 
Eq. ( [2~j~[ >, this corresponds just to replace N c —♦ N c - 1 in the ki¬ 
netic energy contribution. Thus, the projection becomes trivial 
to deal with, and one obtains for a moving baryon of momen¬ 
tum P: 


M b = N c m Q + 


P 2 

2 IIIqNc 


3 (N c - 1 ) + , a2 

2b 2 m Q bx/n 2 


+ b 2 cr), (22) 


where A 2 = zr/12. Minimizing with respect to b (bo) yields 
the baryon mass at rest. At large ihq , bo and the baryon mass 
become: 


bo = 


3 sjn N c — 1 
A 2 m Q N c 

O( l / m o) 


1 -9 n 


N r - 1\ 2 rr A 


N r 


A 6 m 2 

/i iriQj 


M b = N c m Q 


+ (N c - 1 ) 


2 m Q N c 
3 cr I N r 


\ 2 2 niQ 


N r - l 6 n 


2 i4 \ 
Amr 


0(l/m 2 Q ), 


(23) 


which shows a delayed onset of the heavy quark regime due to 
large numerical factors. Thus, one should expect relativity to 
play a role even for moderately heavy quarks. 


3.4. Relativistic corrections 

Of course, a full relativistic treatment implies particle cre¬ 
ation as implied by locality, and Poincare invariant Hamiltonian 
methods with a fixed number of particles exhibit well known 
features (see e.g. Ref. |[28l and included references). While this 
can be improved, here only an estimate of the relativistic cor¬ 
rections is considered by the standard replacement at the single 

particle level, iuq + p 2 /2mQ —» y/r + m 2 Q , which leads re¬ 
markably to an analytical expression for the zero momentum 
projected variational energy 



which reproduces from the simple non-relativistic CM rule 
N c —> N c - 1 in the kinetic energy in the heavy quark limit 0 
The scheme as in the mean field case of minimizing with re¬ 
spect to the oscillator parameter b yields the final baryon mass 
at any N c 0 This case will be used in order to compare with 
the LQCD results in Ref. 00, where the largest quark masses 
used are still not in the heavy regime. 


3.5. Ground state correlations 

As expected Eqs. ((20]) and (19) are N c independent and cor¬ 
respond to the leading order approximation. These equations 
have corrections corresponding to different physical effects. 
Within the Gaussian ansatz for the single particle states Eq. d2T i 


Tn the cr = 0 case one has Mg — 3 m p = -0.00034o^ mg \21\ vs 
Mg ~ 3 mQ = -0.00031 a^mQ from Eq. |2l) . For the case cr ^ 0 more so¬ 
phisticated ansatze were tried embodying better short and long distance behav¬ 
iors, but improvement is at the per cent level since the quarks are located in the 
mid-range region. Discussion of several possibilities will be given elsewhere. 

2 Semiclassical collective quantization methods provide an alternative after 
due attention to zero modes is paid ehh). 


3 Note that here one projects and does not boost the mean field solution. 
In the relativistic case the rest and inertial masses ought to coincide due to 
Poincare invariance. The necessary identity between boosting and projecting 
onto linear momentum only holds for exact solutions G3. At the mean field 
level the identity is guaranteed at the mean fi eld s olution ED. 

4 Note that the direct extrapolation of Eq. |24| to light quarks m q —> 0 leads 
to the rest mass Mb/(N c fd) = 1.81 - 0.50 /N c - 0.19/Af 2 + ..., which is the 
crude estimate for the multiplet center in the quenched approximation. 
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a Harmonic oscillator shell model interpretation applies since 
the baryon is in a (1 s) Nc state. In this picture, ground state cor¬ 
relation correspond to virtual excitations to higher shell states 
(«i/i) ... (n Nc l Nc ). 

In order to quantify the accuracy of the Hartree approxima¬ 
tion within the large N c framework, one evaluates the variance 
of the Hamiltonian defined by A// 2 = (II 2 ) - (H) 2 where 
(O) = (i//\0\i//). When solving the equation approximately, 
as it is done here using a variational wave function, it turns 
out that A H mr /(H) = Oil/ typical of statistical fluctu¬ 
ations. Straightforward calculation, explicitly using the mean 
field equation Eq. (20i, shows that^j 


A Hl = 


N C (N C - 1) 


[< V QQ ') 2 + < Vl Q ,) - 2(V Q qV q , q „)\ . (25) 


Only when the self-consistent Hartree mean field equation is 
exactly satisfied and due to the Casimir scaling assumption, 
Eq. (16 1 , one has A//eartree = 0(n®}, which means A H^/(H) = 
0(1 JN(.) for the correction relative to the baryon mass. 


3.6. Multiquark interactions 

In general, there are multiquark interactions which contribute 
to the baryon mass at the nominal leading O (N c ). For heavy 
quarks one expects that in the baryon only n-body interactions 
with n < N c are of any significance. For N c — 3 there is a 
long history of studying the 3-quark interactions, where there 
are two competing alternatives to confining forces of quarks in 
baryons, the A (pairwise triangle shape) and the Y (junction 
shape) inspired by string models l3T| . 

Three body interactions have been addressed perturbatively 
|32l for arbitrary N c . In the present case, the non-perturbative 
effect of 3-body interactions can be visualized with one exam¬ 
ple. Consider a 3-body potential of the form: 

3 

V 3 (xi , x 2 , X 3 ) = ^ v 3 (x, - X) d abc A a ®A b ® A c , (26) 

;=i 


where X is the CM position of the three quarks. The expectation 
value of V 3 in the baryon ground state at rest can be evaluated 
explicitly choosing v 3 (r) = Aj (-y + cr 3 r) where d 3 and cr 3 

are O (/V'.’j, one obtains for the Gaussian wave function: 


+ +crib )- 


<V 3 > = 2. 

where the color matrix element for the baryon was used, 

4 


(dabc A a ®A b ® A c ) = 4 {Nc , 3)! (M3 ■ 

l\r I 


■5 Nr 


N, 


(27) 


■)■ (28) 


5 Here the notation corresponds to 

(Vqq')eJ d 3 xd 3 y Vqq(x - y)\<p(x)\ 2 \<p(.y)\ 2 
(Vqq' V Q ’Q"^f d 3 xd 3 yd 3 z Vqq(x - y)V QQ (y - z)|^(x)| 2 |^0>)| 2 |^(z)| 2 



m Q 


Figure 2: Baryon mass as a function of the string tension for Gaussian wave 
function. Depicted are the results for non-relativistic (full curves) and relativis¬ 
tic (dashed) calculations, and the lattice QCD results for mo defined by Eq. j2j 
(diamonds) ED The color coding is that of Fig. 1, and in green the limit 
N c —> oo. The string tension corresponding to the lattice QCD results was 
obtained as explained in the text. 


Note that the expectation value of the 2-body interaction 
Eq. p2| ) and the one of the 3-body interaction studied here 
have the same form except that their N c scalings differ by 
terms which are of relative order 1/N 2 . Therefore, the 3-body 
forces cannot be distinguished from the 2-body ones unless 
those higher order terms in the expansion are taken into ac¬ 
count. This is in a sense direct consequence of the mean field 
approximation, which naturally ’’hides” the n-body nature of 
the interactions. Other n-body forces are in principle possi¬ 
ble for a large N c baryon, whose color structure is given by 
l/W"^ 1 d ai ... an A a> ®- ■ -®A a ", where d ai ... aii is the rank7* invariant 
symmetric tensor of S U(N C ). A simple calculation shows that 
they contribute to the baryon mass with an overall factor N c /n\, 
which implies that even for very large N c , 77-body forces with 
77 > 5 become very suppressed. 


3.7. Hype rfine effects 


The simple OGE potential contains hyperfine components 
<9(777g 2 ), which have implications on meson spectra (see e.g. 

Ref. l33l ). as they contribute at O (;V ( (I ) in mesons, but con¬ 
tribute to hyperfine splitting in baryons only at 0(1/N C ). They 
can be easily evaluated as perturbations using the wave func¬ 
tion obtained here. A quick calculation generalizing the N c = 3 
result Il34l to arbitrary N c gives for the hyperfine mass shifts: 


dM H /(S ) - 



a N s c (m Q ) 

m 2 jb 3 


(S(S+1)-^V C ). 


(29) 


They play no role for the spin-weighted average baryon mass 

Eq. Q. 
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Figure 3: Effective potential for a heavy quark in a heavy baryon with iuq = 
2 y/cr (dashed) and 10 y[cr (full). Same color coding as in Fig.pl 


highly interesting to extend the lattice calculations to the non- 
relativistic regime, where the theory can be more easily han¬ 
dled. 

The mean field approximation is visualized through the mean 
field potential V(r) created by the N c - 1 quarks, see Eq. ( fl9| i. 
In the present case, for zero momentum states and the Gaussian 
profile, Eq. (21 1 one obtains: 

b Q cr 


V(r) = 


yf2n 


■ e o 


il 2 -cr 

( h 2 3 

^ + r 2 

- erf 

' V2 r' 

\ 

U J 

r 

{ b 0 J 


(32) 


which is shown for illustration, in Fig. [3] for different values 
of N c and mg. Improvements to this behavior correct for long 
distance behavior and will be discussed in a forthcoming publi¬ 
cation. 


5. Conclusions 


4. Towards relating to LQCD results 

Following the motivation of this work, the aim here is to com¬ 
pare the mean field description including relativistic and CM 
corrections with results from LQCD. At present, the only avail¬ 
able LQCD results for ground state baryon masses at several N c 
values are those of Refs. 0[8j (slightly updated in Ref. Ifl8l ), 
where quenched calculations have been undertaken at several 
values of the quark mass and at N c = 3, 5, and 7. While the 
purpose there was to pursue the light quark limit, here the oppo¬ 
site situation is emphasized where simplifications are expected 
and the quenched approximation is better fulfilled. 

As discussed earlier, the explicit N c dependence is inferred 
from taking cr to be N c independent. The lattice results dis¬ 
played in Refs. 01 m are given in lattice units, with a the 
lattice spacing. Using the form of the quark-quark potential the 
Sommer parameter r\ is determined by the standard definition 

- r\V' QQ {n) = -1 , (30) 

yielding in the present case 

r i cr = i - ^ • (3D 

This value, namely r\ \J7r = 0.859, is roughly valid for the 
LQCD calculations with N c — 3, 5 and 7, where the respec¬ 
tive results from Table I of Ref. 0 are 0.856(5), 0.850(4) and 
0.845(2). Using the values of r\ ju in the same Table one ob¬ 
tains respectively yjcra = 0.219(2), 0.225(2), and 0.216(1). 
For the level of precision of the present comparison it is there¬ 
fore sufficient to take y[cr a = 0.22 for all N c . While the main 
goal of 0® was to pursue the lowest quark mass limit, some 
moderately high quark masses were included. These are now 
used to compare with the results of this work. 

The numerical results are presented in Fig. [2] As expected, 
the relativistic limit sets in at about \[ir ~ m.Q. The lattice 
data of Ref. 0 ® stop at twice larger values, so it would be 


In the present work, a scheme is put forward where the 
large N c expansion of baryon masses in the lattice can be de¬ 
scribed in terms of the mean field approximation as originally 
advocated by Witten and 1 /N c corrections thereof. The quark- 
quark potential is assumed to follow Casimir scaling at arbi¬ 
trary N c and hence proportional to the quark-antiquark poten¬ 
tial, which to good accuracy as per current LQCD calculations 
is /V, -independent. This provides a universal N c independent 
scheme where the ratio of the baryon mass to N c \fir can be 
numerically evaluated. 

It was shown that the corrections to the mean field energy are 
generically Oi yfW c ), but become O(N^), when the mean field 
energy takes its minimum value. This accuracy is the result 
of the density of quarks in the baryon growing as proportional 
with N c . Among the estimated corrections are the leading in N c 
relativistic 0(nr 2 ) and subleading 0(N °) CM corrections. Hy- 
perfine splittings are removed by suitably averaging over spin 
states. When compared with available LQCD calculations, the 
present results account within 20% for the dimensionless ratio 
(Mb ~ N c niQ)/(N c yfcr) which is of natural size. This is encour¬ 
aging, as it suggests to push the LQCD calculations to heav¬ 
ier quark masses and also refine the calculations in the present 
work. 

One of the obvious benefits of the present investigation is the 
possibility of going beyond the ground state and extend these 
ideas to the excited baryon spectrum, where lattice calculations 
are admittedly more involved and less accurate. LQCD calcula¬ 
tions of excited baryons for N c > 3 may still be an unreachable 
goal. However, it is likely that this will be achieved first with 
heavy quarks, and in that case the approach followed here can 
be easily used to predict the excited states. Finally, other heavy 
baryon properties, such as form factors, are easily derived with 
the wave functions obtained here. 
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